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Abstract 

The effects of edge covalent functionalization on the structural, electronic and optical prop- 
erties of elongated armchair graphene nanoflakes (AGNFs) are analyzed in detail for a wide 
range of terminations, within the framework of Hartree-Fock-based semi-empirical methods. 
The chemical features of the functional groups, their distribution and the resulting system 
symmetry are identified as the key factors that determine the modification of structural and op- 
toelectronic features. While the electronic gap is always reduced in presence of substituents, 
functionalization-induced distortions contribute to the observed lowering by about 35-55%. 
This effect is paired with a red shift of the first optical peak, corresponding to about 75% of 
the total optical gap reduction. Further, the functionalization pattern and the specific features 
of the edge-substituent bond are found to influence the strength and the character of the low 
energy excitations. All these effects are discussed for flakes of different width, representing 
the three families of AGNFs. 

Keywords: ZINDO, AMI, substitution, UV-vis spectrum, graphene nanoribbons, configura- 
tion interaction 
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The electronicPfi] optical'^ and magnetic properties'^^ of graphene nanoribbons are closely 
related to their edge morphology and indeed several theoretical proposals have been advanced to 
tune the properties of these systems through edge modifications. For instance, the termination of 
edge C atoms with different groups^^^] and the introduction of magnetic impurities'^^ at the 
edges of zigzag graphene nanoribbons represent promising routes to fully exploit their potential 
for spintronics applications. In armchair graphene nanoribbons (AGNRs), on top of the intrinsic 
tunability given by the well-established family behavior,^ edge modifications can further engineer 
the energy gap and work functionp^^ as well as their optical response.^HUl 

An additional aspect of edge modifications is related to the appearance of structural distor- 
tions: ripples and twists can arise, either localized at the edge region^^or extended to the whole 
structure, with effects on a larger scale.^HSD Indeed, the study of structural modifications has 
drawn itself great attention, since they can lead to a futher tuning of the graphene nanostructure 
properties J 32 * 33 ' We here focus on unraveling the effects of functionalization-induced structural, 
electronic and optical modifications, considering not only the functional group itself but also the 
functionalization pattern. Such analysis can be particularly relevant to design new optoelectronic 
functionalities in graphene nanostructures, above all in light of the recent improvements in their 
controlled synthesis at the nanoscale.^HH 

We carry out our analysis on short armchair graphene ribbons functionalized with different 
edge-substituents, ranging from single atoms to small molecular moieties. Family-dependence is 
also taken into account. We recognize that several factors -including the system symmetry, the 
chemical features of the functional groups and their distribution- are responsible for the appear- 
ance of different structural and optoelectronic modifications. In general, the steric hindrance of 
substituents is found to induce local distortions at the edges, but global effects involving the flake 
backbone can also arise in specific cases. These structural modifications are found to reduce the 
electronic gap up to 300 meV, corresponding to 55% of the total gap decrease. This effect is asso- 
ciated to a similar red- shift of the first optical peak. On the other hand the chemical features of the 
substituents and their arrangement at the edges crucially influence the distribution of the frontier 
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orbitals. Local variations arise as a consequence of the modified molecular symmetry upon edge 
functionalization, while global localization effects are noticed when a large longitudinal dipole 
component is introduced by polar substituents. Moreover, such charge redistribution results in the 
appearance of an additional peak in the low energy region of the spectra, related to the optical 
activation of an otherwise dark excitation. 

Methods and System Description 

The results presented in this paper are obtained within the framework of Hartree-Fock based semi- 
empirical methods,^ which are well tested and reliable for the evaluation of the electronic and 
optical properties of C-conjugated low-dimensional sy stems The AMI modeP is adopted 
for structural optimization (0.4 kcal-mol /A threshold for the forces) and for the study of the 
electronic properties, including the calculation of electric dipole moments. The electronic gap is 
computed as the difference between vertical ionization potential (IP) and electron affinity (EA), 
which are in turn obtained from the total energy of the neutral and (±1) charged states, i.e. IP = 
E(+l) - E(0) and EA = E(0) - E(-l). The optical spectra are evaluated by means of the ZINDO/S 
model,^ with single excitation Configuration Interaction (CIS). Our convergence tests over the 
number of occupied and virtual molecular orbitals (MOs) indicated that a CI energy window of at 
least 4.5 eV below the HOMO and 3.5 eV above the LUMO is required for a reliable characteriza- 
tion of the low energy optical excitations. 

We perform our analysis on elongated graphene nanoflakes (GNFs) with armchair-shaped 
edges, which can be seen as finite portions of graphene nanoribbons. As for infinite ribbons, we in- 
troduce a width parameter N, corresponding to the number of dimeric C lines across the y axis (see 
[l^). According to their width parameter TV, both ribbons^ and flakes 23 present three scaling laws 
for the energy gap, also in presence of edge substituents. This allows us to classify them into three 
different families, following a modulo 3 periodic law: N = 3p + m, with p integer and m = 0, 1 , 2. 
In the following, we will mainly focus on a TV = 7 GNF (width ~ 7.3 A) of length ~ 24 A (x axis, 
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Figure 1 : H-terminated armchair graphene nanoflakes (AGNFs) of fixed length (x axis) and vari- 
able width (y axis), chosen as representative of the three different families of AGNFs, i.e. N=3p, 
3/7+1 and 3p + 2. Edge functionalization is performed by replacing each second H atom along 
the length with a foreign substituent. Two functionalization schemes can be identified: the A-A 
scheme, where the same zigzag line is functionalized on both edges, and the A-B scheme, where 
functionalization involves alternating zigzag lines. Indicated are the point symmetry groups of each 
structure, clean (black), and functionalized in the A-A (red) and A-B (blue) scheme, respectively, 
for ideal (non-distorted and symmetrically terminated) flakes. 



5 



see[jj). We will then extend our analysis to the family dependence by considering two additional 
structures (shown in[l})-c), having widths 8.6 A (N=S) and 9.8 A (N=9) and fixed length. For each 
structure, the lateral ends along the x axis are shaped in such a way to minimize the influence of 
zigzag edges, which are known to affect the electronic properties of the system.E§E2] 

NH 2 COCH 3 



A-A A-B A-A A-B 




Figure 2: Structurally optimized N=7 and N=S AGNFs, functionalizalized with NH 2 (left) and 
COCH 3 (right) groups, according to the A-A and A-B schemes. The angle a is introduced to 
quantify the distortion angle, defined as the maximum deviation from planarity of edge C atoms. 

The effects of edge substituents on the electronic and optical properties of GNFs are investi- 
gated for a number of light-element metal-free functional groups, including atoms and molecular 
groups typically involved in electrophilic aromatic substitutions.^^ We group these terminations 
into three main classes, with respect to the nature of their covalent bond with the graphene edge: 
the first class comprises monovalent halogen terminations (F, CI and Br); the second consists of 
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groups forming a direct bond at the edge, in which either N or O shares alp lone pair with the 
edge C atom (NH 2 , OCH 3 and OCF=CF 2 ); as a third class, we finally consider small moieties 
bonded to the GNF through a C-bond and including a C-double bond, which allows to extend the 
C-conjugation of the graphene backbone for one additional bond length beyond the edge (COCH 3 , 
CH=CH 2 and CH=CF 2 ). The choice of these terminations is basically motivated by their experi- 
mental feasibility in systems similar to those investigated here.'^EIHS] 

Edge covalent functionalization is here performed by replacing each second H atom passivating 
the armchair edges in the reference system, through the regular alternation of substituents. As 
indicated in [TJ two functionalization schemes can be designed: in the A-A scheme the same C 
zigzag line is functionalized on opposite edges, while in the A-B scheme alternating C zigzag 
lines are functionalized. In addition, full substitution, i.e. functionalization of all C edge atoms, is 
considered for monoatomic halogen terminations only, as justified by their reduced size. Note that, 
depending on the chosen scheme and on the flake width, each structure has different symmetry. 
H-terminated flakes have D2h (Cih) symmetry for N odd (even), which is modified by the presence 
of functional groups at the edges. For N odd, the A-A scheme results in a C2 V symmetry, while 
in the A-B scheme the symmetry is C2h- On the contrary, for N even the symmetry is preserved 
as C2h in the A-A scheme and lowered to C s in the A-B scheme. Even though this classification 
is formally correct only for ideal non-distorted and symmetrically terminated structures, since the 
symmetry is always lowered upon structural optimization, we will keep this notation throughout 
the paper for convenience. 

Results and Discussion 

Edge covalent functionalization is known to be an effective strategy to tune the electronic gap 
of graphene nanostructures J ^ * 22 * 32 * As suggested already for simple aromatic compounds,^^ 
the presence of edge functional groups modifies the electronic distribution within the graphene 
network in the edge proximity. As a consequence, heterospecies with different electronegativity 
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than C and covalently bonded at the edges of a graphene ribbon allow to modify the ionization 
potential (IP) and the electron affinity (EA) of the system with respect to the reference H-passivated 
flake. 23 The optical characteristics are then directly affected. 26 However, the edge-decoration of 
the flakes induces also structural modifications, which are closely related and interdependent with 
the electronic stabilization, as we will see in the following. 

[2] shows that, depending on the value of N (even or odd), the type of substituent and the func- 
tionalization scheme (A-B or A- A), different types of distortions can arise, either localized at the 
edges or extended to the flake backbone. For instance, for COCH 3 we have a typical case of 
structural modifications mostly localized at the edges: there is a slight difference for different N 
and scheme, however there are no global, coherent modifications. In this case, local steric effects 
among the functional groups play a primary role in driving distortions: COCH 3 groups, which can 
reorient to minimize the repulsion, induce relatively small tilt of edge C rings to opposite direc- 
tions, perpendicular to the flake basal plane (z axis, see [2]). As a second example, we highlight 
the case of the NH 2 termination: here the interplay of functionalization pattern, group rigidity and 
steric effects can produce coherent distortions involving the whole flake backbone, which thus sta- 
bilize the structure. This is evident for N=7 (N=S) in the A-A (A-B) scheme ([2]), where the system 
tends to twist along its longitudinal axis (x). For N=S in the A-A scheme, a "rippling" mode tends 
also to appear on top of the local distortions induced at the edges. 

To quantify these effects, we introduce an angle a, defined as the maximum deviation of edge 
C atoms with respect to the reference planar configuration (H-GNF, a=0). The values of a for 
different functional groups are reported in[TJfor N=7 GNFs (see Supporting Information for N=S 
and N=9, Table S3). In the case of halogen termination, similar distortion angles are observed 
for both A-A and A-B schemes, while a significantly larger tilting angle is produced upon full 
functionalization. For molecular substituents, where larger distortions are expected in view of 
their size, the tilt of edge C-rings is considerably reduced when the groups are relatively free to 
rearrange, as we described for COCH 3 . This is also the case of methoxy (OCH 3 ) and OCF=CF 2 
groups, where the flexible O-linkage allows them to assume an ordered alignment (see Supporting 
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Information, Figure SI). On the contrary, larger a angles (and global distortions) are noticed 
for rigid groups like NH 2 and ethene units (CH=CF 2 and CH=CH 2 . We remark finally that, 
even when we find discernible structural modifications, they do not result in any significant local 
deviation from the regular graphene bonding and no strongly localized defects appear. For further 
structural details (edge C-C and C-substituent bond lengths), see Supporting Information, Table 
S1-S2. 



Table 1: Distortion angles a (in degrees) for N=7 GNFs functionalized according to the indi- 
cated schemes (I st column). The angle a indicates the maximum deviation of edge C atoms 
as compared to the reference planar H-GNF, as shown in|2| 





Halogen 


Direct bond 


Cbond 




F 


CI 


Br 


NH 2 


OCH 3 


OCF=CF 2 


COCH3 


CH=CF 2 


CH=CH 2 


A-A 


7.2 


11.3 


12.0 


36.5 


11.1 


22.2 


10.5 


19.2 


25.3 


A-B 


6.6 


10.6 


12.5 


22.8 


8.7 


15.3 


12.5 


26.2 


15.5 


full 


17.1 


31.3 


29.1 















We now move to the electronic properties and list in [2] the total dipole arising in presence of 
different substitution patterns. We notice that a large dipole along the flake longitudinal axis (x) 
appears for odd (here N=7) in the A-A functionalization scheme, while it is almost negligible for 
the other scheme. The situation is reversed for ,/V even (e.g. N=S, see Supporting Information, Table 
S4). As shown inland in [3^, the dipole orientation and intensity depend on the characteristics of 
the functional group: electron-donating (e.g. NH 2 ) and -withdrawing groups (e.g. F and COCH 3 ) 
produce large longitudinal dipoles oriented in opposite directions (see also [3^). The presence of a 
sizeable longitudinal dipole (A-A pattern for A^=7 GNF) is accompanied by a spatial localization of 
the frontier orbitals on opposite sides of the flake. As depicted in [3]), for jlx^O the potential well 
describing the flake assumes a "saw-tooth" profile (see e.g. Ref.^), whose height and orientation 
depends on the magnitude and on the sign of \x x . Hence the HOMO and LUMO states present a 
similar asymmetric character in both cases, but with spatially opposite localization (see [3})). On 
the other hand, for the A-B scheme the % distribution of the frontier orbitals undergoes only local 
modifications, specifically related to the altered symmetry of the flake,^2l as shown in |3p. The 
modified symmetry of the frontier orbitals is particularly evident for the NH 2 termination, while it 
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is less apparent for groups bonded through a C bond, such as COCH 3 . See Supporting Information 
for a more complete analysis of the substituents and family dependence (Figure S1-S2). 

Table 2: Components of the dipole moment for edge-functionalized N=7 GNFs (one sub- 
stituent per class), according to the A-A and A-B functionalization schemes. 



Group 


Scheme 


Hx [D] 






F 


A-A 


4.848 


0.479 


0.028 


A-B 


-0.001 


-0.003 


-0.061 


NH 2 


A-A 


-9.335 


-0.026 


-0.121 


A-B 


0.339 


-1.238 


3.823 


COCH3 


A-A 


9.764 


1.138 


-0.049 


A-B 


0.166 


-0.139 


0.269 




Figure 3: (a) A-A functionalization of N=7 GNF by means of electron-donating (e.g NH 2 ) and 
-withdrawing (e.g. COCH 3 ) groups gives rise to longitudinal dipoles (fi x ) oriented in opposite 
directions. The frontier orbitals of NH 2 - and COCH 3 -substituted N=7 GNFs are shown for the (b) 
A-A and (c) A-B schemes. The square well potential profiles modulating the molecular orbital 
distribution are sketched in presence (b) and absence (c) of a longitudinal dipole along the flake. 

In [3] we report the differences of IP and EA (AEA and AIP) of N=7 functionalized GNFs with 
respect to the reference hydrogenated flake, as well as their energy gap (Eq). As expected, an 
uneven shift of IP and EA is produced upon functionalization, related to the electronegativity of 
the substituents; for instance, highly electronegative halogen terminations upwards shift EA more 
significantly than IP, while strongly electron-donating amino groups downwards move IP more 
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Table 3: Differences of EA and IP (AEA and AIP, in eV) with respect to the reference hy- 
drogenated (H) flake and energy gap (Eg, in eV) computed through AMI model for N=7 
GNFs, functionalized according to the A-B scheme. For selected groups, we report the values 
of AEA, AIP and Eg also for the A-A scheme and for H-terminated model flakes with the 
corresponding distorted geometries (in italics). 





H 


Halogens 


Direct bond 


C-bond 


F 

A-A A-B 


CI 


Br 


NH 2 
A-A A-B 


OCH 3 


OCF=CF 2 


COCH 3 
A-A A-B 


CH=CF 2 


CH=CH 2 


AEA 


0.00 


0.57 0.58 
0.01 0.01 


0.53 


0.62 


-0.18 -0.14 

0.08 0.04 


-0.12 


0.70 


0.74 0.69 
0.01 0.01 


0.67 


0.14 


AIP 


0.00 


0.37 0.32 
-0.10 -0.10 


0.30 


0.39 


-0.75 -0.79 
-0.24 -0.19 


-0.50 


0.31 


0.49 0.42 
-0.10 -0.10 


0.19 


-0.22 


E G 


5.13 


4.92 4.87 
5.01 5.01 


4.90 


4.90 


4.55 4.47 
4.81 4.90 


4.75 


4.74 


4.88 4.86 
5.01 5.01 


4.65 


4.76 



than EA. In any case, a gap reduction is always encountered. In the case of C-bonded terminations, 
a further contribution to the gap reduction is given by the increased effective width of the flake, 
driven by the extended ^-conjugation. 

In addition to this well-established behavior, we here more directly estimate the role played 
by structural distortions. To quantify this effect, we consider model flakes obtained by remov- 
ing the functionalizing groups from the final structure, and by then H-passivating the remaining 
dangling bonds: this model structure allows us to decouple the influence of distortions from the 
overall influence of edge functionalization. By inspecting [3j we find a total gap reduction upon 
functionalization ranging from 200 to 700 meV. The comparison with the values of Eq obtained 
for the model distorted H-GNFs (in italics in[3]) points to a non-negligible contribution of structural 
distortions, corresponding to 35-55% of the total amount. These results are very similar to what 
expected for thermal deviations from the ideal structure. 60 It is also worth noting that this trend 
is mostly independent of the functionalization scheme, with a maximum difference of 80 meV for 
the NH 2 case, where the two schemes induce significantly different distortion patterns. 

Starting from the AMI results described above, we now turn to the analysis of optical exci- 
tations, as obtained from ZINDO calculations. To address the role of both molecular-symmetry 
and structural modifications, we first consider the N=l GNF functionalized with NH 2 , where these 
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0.0 2.0 2.5 3.0 3.5 

Energy [eV] 



Figure 4: Optical spectra of GNFs of width parameter N=l , functionalized with NH 2 groups, ac- 
cording to the (a) A- A and (b) A-B functionalization schemes. In addition to the spectrum of the 
planar H-terminated flake taken as reference (grey shaded area, with the main excitations high- 
lighted in a), we show the spectrum of model GNFs obtained by removing the edge substituents 
from the final structures and H-passivating the dangling bonds (dashed lines), in order to perform 
a rigorous analysis of distortions effects. All the curves are obtained through a Lorentzian broad- 
ening of 50 meV. 
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effects are both prominent. In [4] we report the optical spectra for both functionalization schemes, 
A-A (red) and A-B (blue), together with the spectrum of the flat H-GNF (grey shaded area), taken 
as a reference. In addition, we include the spectrum of the H-terminated distorted model flake 
described above (dashed line), which enables a direct analysis of the distortion effects. As already 
discussed,^ the optical spectrum of the planar N=7 H-GNF shows two main peaks in the low- 
energy region with longitudinal polarization (along the x axis), labelled as PI and P2; in addition, 
a transversally polarized weak excitation (PO) is found at lower energy. Our results indicate that 
edge functionalization is responsible for a red shift of the PI peak by 380 and 510 meV in the A-A 
(|4^) and A-B schemes Q)), respectively. Remarkably, a further peak is observed in the optical 
spectrum of the A-B functionalized flake, related to the oscillator strength (OS) increase of PO by 
about 3 orders of magnitude. This OS gain is a purely symmetry effect: the molecular symmetry 
reduction induced by functionalization (from D2/, in the H-GNF to in the A-B scheme) implies 
that the HOMO state acquires the same symmetry as the HOMO-1 one, so that the HOMO — > 
LUMO transition is now allowed to contribute to the PO excitation (see [4]). This symmetry effect 
is in fact not observed in the A-A functionalized flake. 

Concerning the role of structural distortion, we observe that the largely distorted geometry pro- 
duced by the A-A functionalizion (dashed line in [4^) is responsible for about 75% of the optical 
gap reduction, as noticed by comparing the excitation energy of PI for the three curves in [4^. PI is 
now energetically degenerate with the dark excitation PO, also red shifted of about 200 meV with 
respect to the reference H-GNF (see [4]). It is worth noting that the frontier orbital localization ob- 
served for the A-A scheme in presence of polar functional groups (see|3j5) does not affect the basic 
features of the main excitations: despite their spatial localization, the overlap between the frontier 
orbitals is still large enough to preserve the intensity of PI, in comparison with both the reference 
H-terminated and the A-B functionalized GNF In the case of the A-B functionalization pattern, 
distortions contribute to about 45% of the whole optical gap red shift: the impact of distortions is 
less pronounced for this structure, and molecular symmetry effects dominate (see[4}5). 

We next inspect the role of the chemical specificities of edge substituents, by comparing the 
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Table 4: Energy, oscillator strength (OS) and composition in terms of molecular orbital tran- 
sitions of the main excitations of functionalized N=l GNFs, in presence of one prototypical 
substituent per class, i.e. F, NH 2 (also A-A scheme) and COCH 3 , in addition to the reference 
hydrogenated system. In the last column we report the MO transitions contributing to each 
excitation which have a relative weight larger than 0.1 



Termination 


Excitation 


Energy [eV] 


OS 


Transitions (weight) 


H 


P0 


2.54 


0.0005 


H-l -»L(0.37) 
H ->■ L+l (0.39) 


P1 


2 72 


1 62 


H — 5> T (0 79 s ) 


P2 


3.32 


5.43 


H-3 ^L+3 (0.12) 
H-l ->■ L+l (0 76) 


F 


P0 


2.43 


0.21 


H-l -»L (0.27) 
H-»L(0.12) 
H ->■ L+l (0.37) 


PI 


2.58 


1.54 


H ->■ L (0.68) 


P2 


3.26 


5.11 


H-3 ->■ L+3 (0.20) 
H-l ^L+l (0.71) 


NH 2 - A-A 


PI 


2.34 


1.98 


H ->■ L (0.75) 


P0 


2.35 


0.02 


H-2 ->■ L (0.29) 
H ->■ L+l (0.39) 


P2 


3.24 


3.59 


H-2 ->■ L+l (0.50) 


NH 2 - A-B 


PI 


2.21 


1.29 


H ->■ L (0.63) 


P0 


2.48 


0.70 


H-2 ->■ L (0.22) 
H-»L(0.18) 
H ^ L+l (0.40) 


P2 


3.22 


3.92 


H-2 ->■ L+l (0.53) 
H^L+1 (0.10) 


COCH3 


P0 


2.42 


0.04 


H-l ^L(0.34) 
H ^ L+l (0.40) 


PI 


2.55 


1.68 


H -> L (0.78) 


P2 


3.21 


5.13 


H-l -> L+l (0.71) 
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0,0 2.0 2.5 3.0 3.5 

Energy [eV] 



Figure 5: Optical spectra N=7 GNFs functionalized in the A-B scheme with each terminations 
considered in this work, including (a) halogen atoms, (b) directly- and (c) C-bonded groups. In 
each panel, we indicate for one of the groups the three main excitations described in the text, i.e. 
P0, PI and P2. Moreover, the spectrum of the reference H-flake is reported (grey shaded area in 
each panel). All the curves are obtained through a Lorentzian broadening of 50 meV. 



15 



optical spectra of N=7 GNFs functionalized with all the anchoring groups considered in this work, 
as shown in [5] We here consider only the A-B pattern, since negligible differences are observed 
between the two schemes for most of the substituents (see Figure S3 in the Supporting Information 
for more details). The analysis of fully halogen terminated GNFs is also included in the Support- 
ing Information (see Figure S4 and Table S8). A non rigid red shift of the spectra is generally 
observed upon functionalization, with a larger red shift for PI than for P2. All the spectra of the 
groups anchored through a direct bond ([5J5) display the same key feature already discussed for the 
amino termination and reported in[4j): an additional distinct peak emerges, corresponding to the 
PO excitation which becomes optically active. This effect is also present, even though less intense, 
for halogen terminations: PO appears as a shoulder in the main PI peak for fluorinated flakes (|5^). 
In the case of C-bonded molecular substituents, the impact of functionalization on the electronic 
states is further reduced, so that the main spectral features are basically preserved with respect to 
the flat H-terminated flake (see spectra in [5}) and the composition of the main excitations in [4]). 



(a) N=8 (b) N=9 




Energy [eV] Energy [eV] 

Figure 6: Optical spectra of symmetrically functionalized GNFs of width parameter (a) N=S (A-A 
scheme) and (b) N=9 (A-B scheme). One edge substituent for each class of terminations considered 
in this work is addressed for both GNFs. The spectrum of the reference H-GNF is reported in the 
background of each panel (grey shaded area) with the indication of the main excitations (PO, PI 
and P2). All the curves are obtained through a Lorentzian broadening of 50 meV. 

In order to complete our analysis, we finally investigate the spectra of the GNFs of width pa- 
rameter N = & and N = 9, functionalized with one prototypical group per class, namely F, NH 2 
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and COCH3, according to the higher symmetry scheme (A- A for N=S and A-B for N=9). The 
optical spectra computed for these systems are shown in|6j with the reference spectrum of the cor- 
responding planar H-GNF reported in the background of each panel (grey shaded area). The optical 
spectrum of the iV=8 GNF is dominated by a first intense peak PI, which is about 1 eV lower in 
energy compared to those of the other GNF families (compare [6^ with [5] and |6j)), in agreement 
with the trend for the electronic gaps. 23 Contrary to N =7, this is always the lowest energy exci- 
tation, regardless the edge termination (See Supporting Information, Table S5). Moreover, we do 
not observe any functionalization-induced modification in the it distribution of frontier orbitals, 
since the molecular symmetry remains the same in the A-A scheme (Q-2\x)- In this case the main 
effect of functionalization is thus to red shift the spectrum with respect to the reference H-GNF, in 
an almost rigid fashion (up to over 200 meV upon amino termination). 

On the other hand, the optical spectra of N=9 GNFs manifest more similarities with those 
of N=7 GNFs (see [6Jd and [5] for comparison), both in the spectral shape and in the composition 
of the main excitations (see Supporting Information, Table S6). Also for the ribbons belonging 
to the N=3p family, the P0 excitations becomes optically active, giving rise to a distinct peak in 
presence of NH 2 functional groups, according to the mechanism described for the N=l GNFs. 
The lowest energy peak is red shifted by about 300 meV upon amino substitution, compared to 
the reference H-GNF. In presence of other edge substituents (F and COCH 3 ), the energy of PI is 
basically unaffected, while P2 is red shifted of about 200 meV upon each termination (see [6^). 
This is an effect of the MO distribution observed for the GNFs belonging to this family, for which 
the frontier orbitals present reversed orientation of the % distribution with respect to the flakes of 
the N=3p+l family (see Figure S2). 

Conclusions 

We have analyzed in details the effects of edge covalent functionalization on the structural, elec- 
tronic and optical properties of armchair graphene nanoflakes. Different functionalization groups 
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and patterns have been discussed for a specific flake, and the analysis extended and rationalized 
for additional systems, according to the armchair width modulation. Our results indicate that 
steric repulsion among substituents generally determines local distortion at the edges, even though 
structural modifications can involve the whole flake backbone in specific cases. Functionalization- 
induced distortions are found to play an effective role in the electronic gap reduction experienced 
by the flakes, with contributions ranging from 35 to 55% of the total gap decrease. Depending on 
the functionalization pattern and the character of the substituents, interesting charge-redistribution 
effects are seen for the frontier orbitals, either related to the modified molecular symmetry or to the 
appearance of a large longitudinal dipole component. All these effects are reflected in the optical 
spectra: distortions largely contribute to the lowering of the optical gap (up to 300 meV), while the 
appearance of an additional peak at low energy is related to the distribution and the characteristics 
of the functional groups. Our results indicate that the interplay among the system symmetry, the 
chemical features and the distribution of the substituents, as well as the resulting structural distor- 
tions is crucial to understand the modification of the opto-electronic properties of functionalized 
graphene nanoflakes. 
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The Supporting Information is organized in three sections. In the first one we report the struc- 
tural details (bond lengths and distortion angles) of the optimized N=l , N=S and N=9 graphene 
nanoflakes (GNFs), to complete the data reported in the main text. The second section supports 
the analysis on the electronic properties, including a table with the dipole moment components of 
N=S and N=9 GNFs (functionalized according to both A-A and A-B schemes) and the isosurfaces 
of the frontier orbitals for the considered systems. The last part is dedicated to optical properties. 
The analysis on the effects related to different functionalization schemes and structural distortions 
are extended here also to COCH 3 and F terminations, as representatives of the classes of C-bonded 
and halogen substituents. We also include the tables with the composition of the main excitations 
for N=S and N=9 GNFs, and the analysis of the optical features of full halogen-terminated GNFs 
(only N=T). 



Structural Properties 

We report in Table SI and Table S2 the bond lengths - C-C and C-X bonds, respectively, with 
X=F, CI, Br, N, O, C, depending on the functional group - at the edge of functionalized graphene 
nanoflakes (GNFs). 

o 

Table SI: C-C bond length (in A) at the edge of each considered functionalized GNF, opti- 
mized with AMI model. The case of H termination is also reported as a reference. 



C-C bond length 



H 


Halogens 


Direct bond 


C-bond 


F 


CI 


Br 


NH 2 


OCH 3 


OCF=CF 2 


COCH3 


CH=CH 2 


CH=CF 2 


1.40 


1.38 


1.38 


1.38 


1.40 


1.38 


1.37 


1.38 


1.38 


1.38 



Table S2: Bond lengths (in A) between edge C atoms and functional groups for each consid- 



Halogen 


Direct bond 


Cbond 


F 


CI 


Br 


N 





C 


OCH3 


OCF=CF 2 


COCH3 


CH=CH 2 


CH=CF 2 


1.39 


1.70 


1.88 


1.40 


1.38 


1.40 


1.49 


1.46 


1.45 



In Table S3 we report the values of the distortion angle a for N=S and N=9 functionalized 
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GNFs with selected substituents (F, NH 2 and COCH 3 ) for both A-A and A-B schemes. As intro- 
duced in the main text, a is defined as the maximum deviation from planarity of edge C atoms 
in functionalized GNFs. By inspecting Table S3 we notice that in general, as discussed for N=l 
GNFs in the main text, larger distortions are noticed also for N=S and N=9 GNFs in the schemes 
corresponding to lower symmetry models (A-B and A-A, respectively). 

Table S3: Distortion angle a, expressed in degrees, of optimized GNFs of width parameter 
Af=8 and Af=9, functionalized with selected groups (F, NH 2 and COCH) according to both 
schemes. For F termination, the case of full functionalization is also addressed. The angle 
a represents the maximum deviation of edge C atoms from the planar configuration of the 
reference H-termi nated flake. 



A-A 
A-B 



N=9 



7.2 
36.8 
~9XT 



NH 2 
~T7~8~ 
38.4 



COCH, 
12.5 



99 9 



9.4 
4.9 
~2L4~ 



NH 2 
24.7 
10.1 



COCH, 
9.3 " 
13.5 



Electronic Properties 

In Table S4 we report the components of the dipole moments of functionalized GNFs of width 
parameters N=S and N=9 for both A-A and A-B scheme. As discussed in the main text, we notice 
that in general a non negligible dipole component along the longitudinal axis (/i x ) appears in the 
scheme corresponding to the lower nominal symmetry (A-B for /V=8 and A-A for N=9). The 
sign of \l x is determined by the electronegativity of the substituents. The results are basically in 
agreement with those discussed in the main text for N=7. 
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Table S4: Dipole components of funcitonalized GNFs of width parameter N=8 and N=9: one 





Group 


ocneme 


11 i n 1 


11 i n 1 


11 i n 1 
I4->J 




F 


A-A 


005 


-0 001 


-0 81 Q 




A-R 


1 91? 


9QQ 


675 


oo 

II 


NH 2 


A-A 


-0 960 


-0 081 


116 




A-R 


-16 317 


-0 168 

V7. 1VJO 


-0 6Q1 




COCH3 


A-A 


-1 Q64 


1 516 


-0 414 




A-B 


8.190 


0.240 


-0.191 




F 


A-A 


7.409 


0.002 


0.076 




A-B 


0.000 


0.008 


0.187 


II 


NH 2 


A-A 


-2.815 


0.260 


-2.321 




A-B 


-1.191 


0.970 


1.594 




COCH3 


A-A 


9.806 


0.034 


0.200 




A-B 


-2.128 


-1.914 


0.105 



We present the isosurfaces of the frontier orbitals of N=l (Figure SI), A^=8 and Af=9 GNFs 
(Figure S2) for the functionalized systems analyzed in the main text. For consistency, A-B (A-A) 
scheme is addressed for flakes with odd (even) N. 
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LUMO 
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Br 





CH=CF 



CH=CH 



Figure SI: Isosurfaces of the frontier orbitals of functionalized GNFs (N=l) according to the A-B 
scheme. The orbitals of the reference H-terminated system are included for comparison. 
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N=8 

H F NH 2 COCHj 



LUMO 



HOMO 




N=9 

H F NH 2 COCH 3 




Figure S2: Isosurfaces of frontier molecular orbitals of functionalized GNFs of width parame- 
ters N=S (top) and Af=9 (bottom). The frontier orbitals of the reference H-terminated system are 
included for comparison. The A-A (A-B) scheme is addressed for N=S (N=9). 
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Optical Properties 

To supplement the analysis on functionalization-induced distortion effects on the optical properties 
of GNFs, we include in Figure S3 the optical spectra of N=l GNFs, functionalized by means of 
COCH3 groups [panels (a-b)] and F atoms [panels (c-d)]. In addition to the spectra of functional- 
ized GNFs (solid lines), we include also those of the H-terminated distorted model structures, ob- 
tained by removing edge substituents from the final structure and then H-passivating the remaining 
dangling bonds (dashed lines); the spectrum of the reference H-terminated GNF is represented in 
the background of each panel (grey shaded area). Unlike NH 2 -functionalization (see also Figure 4 
in the main text), here the spectra corresponding to the two schemes are basically identical for both 
substituent species, as suggested also from the similar geometry obtained in the two cases (see 
main text, Figure 2b for COCH 3 , and the values of the distortion angle a in Table 1). Compared 
to the reference H-GNF, we notice for all the spectra a red shift of PI (about 150 meV) and of P2 
(about 100 meV). Most part of these red shift effects (about 80%) is to be ascribed to structural 
distortions, as shown by the dashed spectra of the H-terminated distorted model flakes. 
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Energy [eV] Energy [eV] 



Figure S3: Optical spectra of GNFs of width parameter N=7, functionalized with COCH 3 (a-b) 
and F (c-d), according to the A-A and A-B schemes. In addition to the spectrum of the planar 
H-terminated flake taken as reference [grey shaded area, with the main excitations highlighted in 
(a)], we show the spectrum of model GNFs obtained by removing the edge substituents from the 
final structures and H-passivating the dangling bonds (dashed lines), in order to perform a rigorous 
analysis of distortions effects. 
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In Table S5 and Table S6 we report the composition of the main excitations of functionalized 
GNFs of width parameters N=S and N=9, respectively. The corresponding spectra are shown in 
the main text (Figure 6). 

Table S5: Energy, oscillator strength (OS) and composition of the two lowest energy exci- 
tations of functionalized GNFs of width parameter N=8 (A -A scheme, see spectra in Figure 
6a of the main text). In the last column only the molecular orbital transitions with relative 
weights larg er than 0.1 are reported. 



System 


Excitation 


Energy [eV] 


OS 


Transitions (weight) 


H 


PI 


1.84 


3.33 


H->L(0.91) 


P2 


3.64 


3.47 


H-2 ->■ L+2 (0.56) 


F 


PI 


1.76 


3.22 


H->L(0.91) 


P2 


3.56 


3.62 


H-2 ->■ L+2 (0.62) 


NH 2 


PI 


1.69 


3.04 


H-»L(0.91) 


P2 


3.04 


1.64 


H-2 ->■ L (0.46) 
H^L+2 (0.18) 
H^L+4(0.15) 


COCH 3 


PI 


1.73 


3.22 


H-»L(0.91) 


P2 


3.14 


1.08 


H-2 ->■ L+l (0.23) 
H-2 -»L (0.10) 
H ->• L+4 (0.26) 
H^L+6 (0.13) 
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Table S6: Energy, oscillator strength (OS) and composition of the two lowest energy exci- 
tations of functionalized GNFs of width parameter N=9 (A-B scheme, see spectra in Figure 
6b of the main text). In the last column only the molecular orbital transitions with relative 
weights larg er than 0.1 are reported. 



System 


Excitation 


Energy [eV] 


OS 


Transitions (weight) 




DA 

r(J 


/in 
Z.4L) 


A AA 
U.UU 


H-l ^L(0.37) 
H -> L+l (0.41) 


H 


PI 


2.47 


1.71 


H^L(0.85) 




P2 


3.24 


4.47 


H-2 ^L+2 (0.12) 
H-l L+l (0.73) 


F 


P0 


2.34 


0.06 


H-l ^L(0.36) 
H ->■ L+l (0.39) 


PI 


2.48 


1.31 


H ->■ L (0.78) 




P2 


3.12 


4.80 


H-l L+l (0.69) 




PI 


2.17 


0.76 


H ->■ L (0.55) 
H^L+1 (0.12) 


NH 2 


P0 


2.44 


0.50 


H-l ^L(0.27) 
H ->■ L (0.24) 
H ->■ L+l (0.33) 




P2 


3.03 


4.64 


H-l L+l (0.57) 


COCH3 


P0 


2.33 


0.01 


H-l ^L(0.37) 
H ->■ L+l (0.41) 


PI 


2.48 


1.21 


H ->■ L (0.80) 




P2 


3.08 


4.55 


H-l -> L+l (0.64) 
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We finally present the results on the optical properties of fully functionalized N=7 GNFs with 
halogen terminations. In Figure S4(a) we show the computed spectra, where the main excitations 
are indicated for the reference H-GNF (grey shaded area). In case of full functionalization the 
model symmetry of the H-terminated flake is preserved and hence the character of the frontier 
orbitals [see Figure S4(b)] as well as the composition of the main excitations in terms of molecular 
orbital transitions and their intensity are not significantly affected (see Table S7). A red shift of the 
spectra (about 350 meV for PI and about 250 meV for P2) is noticed, compared to the reference 
H-GNF. 



( a ) ( b ) HOMO LUMO 




Energy [eV] 



Figure S4: (a) Optical spectra of N=7 graphene nano flakes fully functionalized with halogen atoms. 
For comparison, the spectrum of the reference H-terminated flake is included (grey shaded area 
with the labels of the main excitations), (b) Frontier orbitals of the corresponding structures. 
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Table S7: Energy, oscillator strength (OS) and composition of the main excitations (PO, PI 
and P2) of fully functionalized N=7 GNFs with halogen terminations (F, CI, Br - see Figure 
S5). In the last column only the molecular orbital transitions with relative weights larger 
than 0.1 ar e reported. 



System 


Excitation 


Energy [eV] 


OS 


Transitions (weight) 


H 


PO 


2.54 


0.0005 


H-l ^L(0.37) 
H — >■ L+l (0.39) 


PI 


2.72 


1.62 


H ->• L (0.79) 


P2 


3.32 


5.43 


H-3 ^L+3 (0.12) 
H-l L+l (0.76) 


F 


PO 


2.36 


0.02 


H-l ^L(0.30) 
H ^ L+l (0.45) 


PI 


2.45 


1.81 


H -> L (0.82) 


P2 


3.21 


4.70 


H-3^L+3 (0.11) 
H-2 -+ L+l (0.73) 


CI 


PO 


2.27 


0.01 


H-l -+L (0.28) 
H -+ L+l (0.42) 


PI 


2.38 


1.62 


H-»L(0.81) 


P2 


3.08 


4.42 


H-3 ^L+2 (0.10) 
H-l -+ L+l (0.67) 


Br 


PO 


2.27 


0.01 


H-l -+ L (0.29) 
H ^ L+l (0.42) 


PI 


2.39 


1.59 


H-»L(0.81) 


P2 


3.06 


4.33 


H-3 ^L+2 (0.10) 
H-l ^ L+l (0.69) 
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